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ABSTRACT 

We investigate the dependence of physical properties of galaxies brighter than Mr = — 18.0+51og/i in 
the Sloan Digital Sky Survey (SDSS) on environment, as measured by local density using an adaptive 
smoothing kernel. We find that variations of galaxy properties with environment are almost entirely 
due to the dependence of morphology and luminosity on environment. Because galaxy properties 
depend not only on luminosity but also on morphology, it is clear that galaxy properties cannot be 
determined solely by dark halo mass. When morphology and luminosity are fixed, other physical 
properties, such as color, color-gradient, concentration, size, velocity dispersion, and star formation 
rate, are nearly independent of local density, without any break or feature. The only feature is 
the sharp decrease of the late type fraction above the critical luminosity of about Mr = —21.3 in 
the morphology versus luminosity relation. Weak residual dependences on environment include that 
of the color of late types (bluer at lower density) and of the L-a relation of early types (larger 
dispersion at higher density for bright galaxies). The fraction of galaxies with early morphological 
type is a monotonically increasing function of local density and luminosity. The morphology-density- 
luminosity relation, as measured in this work, should be a key constraint on galaxy formation models. 
We demonstrate that the dependence on environment of the morphology of galaxies originates from 
variations in density on effective Gaussian smoothing scales much smaller than 12ft,~^Mpc. We find 
that galaxy morphology varies both with density measured on an effective Gaussian smoothing scale 
of 4.7ft.~^Mpc and with distance to the nearest bright galaxy, particularly when the distance is about 
0.2 Mpc. We propose as a mechanism that the morphology of galaxies in galaxy systems is 
transformed by the tidal force. 

Subject headings: galaxies:clusters:general - galaxies:evolution - galaxies:formation - galax- 
ies:fundamental parameters - galaxies: general - galaxies:statistics 



L INTRODUCTION 

The origin of the morphology of galaxies is one of the 
mysteries of galaxy formation. Recently, there has been a 
great deal of effort to understand the origin of morphol- 
ogy by inspecting dependences of galaxy properties on 
environment. This line of study dates back to the 1930's 
(e.g. Hubble & Humason 1931) when it was realized 
that clusters were dominated by ellipticals and lenticu- 
lars and that environmental factors played an important 
role in determining the morphology of galaxies. Oem- 
ler (1974) found the morphology-radius relation; the late 
type galaxy fraction decreases with radius within a clus- 
ter. This relation was confirmed by Dressier (1980) who 
argued that the fraction of morphological types is a func- 
tion of local galaxy density. Postman & Geller (1984) 
extended this morphology-density relation down to the 
group environment. Gunn & Gott (1972) argued that SO 
galaxies in great relaxed clusters like Coma were the re- 
sult of spirals being stripped of gas by ram pressure strip- 
ping due to hot intra-cluster gas. Gott & Thuan (1976) 

^ Korea Institute for Advanced Study, Dongdaemun-gu, Seoul 
130-722, Korea 

^ Department of Physics, Drexel University, 3141 Chestnut 
Street, Philadelphia, PA 19104, USA 

^ Department of Astrophysical Sciences, Peyton Hall, Princeton 
University Princeton, NJ 08544-1001, USA 

* Center for Cosmology and Particle Physics, Department of 
Physics, New York University, 4 Washington Place, New York, NY 
10003, USA 



argued that elliptical galaxies were produced by larger 
initial density fluctuations which would have higher den- 
sity at turn-around and where star formation would be 
completed before collapse. Such larger initial density 
fluctuations would be more likely in a region that would 
later turn into a high density environment. 

Environmental dependences of galaxy properties other 
than morphology have also been studied extensively. Re- 
cently, as large galaxy redshift surveys like the Two De- 
gree Field Galaxy Redshift Survey (2dFGRS; CoUess et 
al. 2001) and the Sloan Digital Sky Survey (SDSS) have 
been completed, extensive studies have been made to ac- 
curately measure the environmental effects on various 
physical properties of galaxies (Goto et al. 2003; Balogh 
et al. 2004a,b; Rojas et al. 2004, 2005; Hoyle et al. 2005; 
Blanton et al. 2005a; Croton et al. 2005; Tanaka et al. 
2005; Weinmann et al. 2006; Boselh & Gavazzi 2006; 
among many others). Works done by others on simi- 
lar topics will be discussed in each section below when 
the dependence of each of physical properties of galaxies 
on environment is studied. The observed environment- 
galaxy property relations provide only indirect clues to 
the origin of morphology, but they certainly are impor- 
tant constraints on models for the physics of galaxy for- 
mation. 

In a companion paper we study the relations among 
various physical properties of galaxies in the SDSS spec- 
troscopic sample (Choi et al. 2006, hereafter Paper I). 
In this paper we extend our work to the environmcn- 
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tal dependence of the physical parameters. We accu- 
rately measure the morphology-dcnsity-lurninosity rela- 
tion from a series of volume-limited samples drawn from 
the SDSS. We also use morphological subsets of early and 
late type galaxies to measure the dependences of various 
physical properties of galaxies on environment once the 
morphology and luminosity arc fixed. By studying the 
environmental dependence at different smoothing scales, 
we address the important question of whether galaxy 
morphology depends primarily on the large-scale envi- 
ronment in which the galaxy initially formed or on a 
smaller scale environment that may reflect the influence 
of later evolutionary effects such as galaxy-galaxy inter- 
actions. 

2. OBSERVATIONAL DATA SET 
2.1. Sloan Digital Sky Survey Sample 

We use a large-scale structure sample, DR4plus, of the 
SDSS (York et al. 2000; Blanton et al. 2003a; Fukugita 
et al. 1996; Gunn et al. 1998, 2006; Hogg et al. 2001; 
Ivezic et al. 2004; Lupton et al. 2001; Pier et al. 2003; 
Smith et al. 2002; Stoughton et al. 2002; Tucker et 
al. 2006) from the New York University Value-Added 
Galaxy Catalog (NYU-VAGC; Blanton et al. 2005b). 
This sample is a subset of the recent SDSS Data Release 
5 (Adelman-McCarthy et al. 2007). Our major sample of 
galaxies used here is a subset of the large-scale structure 
DR4plus (LSS-DR4plus) sample referred to as "voidO," 
which includes the Main galaxies (Strauss et al. 2002) 
with apparent magnitudes in the range 14.5 For the ap- 
proximately 6% of targeted galaxies that lack a measured 
redshift because of fiber collisions, we assign the redshift 
of the nearest neighbor. 

Completeness of the SDSS is poor for bright galaxies 
with r < 14.5 because of spectroscopic selection crite- 
ria (which exclude objects with very large flux within 
the fiber aperture) and difficulties of automatically mea- 
suring photometric properties of very extended sources. 
Due to the magnitude limits of the voidO sample, the full 
range of magnitude of the sample is only AM = 3.1 at a 
given redshift, and is even smaller than that in the case 
of a volume-limited sample with a finite range in red- 
shift. As described in detail in section 2.1 of Paper I, we 
add the missing bright galaxies and thereby extend the 
magnitude range by using various existing redshift cat- 
alogs which we match to the SDSS data. In total, 5195 
bright galaxies are added to the voidO samples, which 
means our whole sample includes 317533 galaxies within 
our angular sample boundaries (see Fig. 1). Volume- 
limited samples derived from the resulting catalog have 
nearly constant comoving number density of galaxies in 
the radial direction for redshifts z > 0.025. We treat 
our final sample as effectively having no bright limit at 
redshifts greater than z = 0.025. More details about this 
sample can be found in Paper I. 

2.2. Volume-limited Samples 

To study the effects of environment on galaxy prop- 
erties it is advantageous for the observational sample to 
have the lowest possible surface-to-volume ratio, to min- 
imize the boundary effects (see Park et al. [2005] re- 
garding similar effects in topology analysis). For this 
reason we trim the DR4plus sample as shown in Figure 



1, in which gray lines delineate our sample boundaries 
in the survey coordinate plane [(A, 77); Stoughton et al. 
2002]. There are also three stripes in the Southern Galac- 
tic Cap observed by SDSS that we do not use because of 
their narrow angular extent. These cuts leave approxi- 
mately 4464 square deg^ in the survey region. Within 
our sample boundaries, we account for angular varia- 
tion of the survey completeness by using the angular se- 
lection function defined in terms of spherical polygons 
(Hamilton & Tegmark 2004), which takes into account 
the incompleteness due to mechanical spectrograph con- 
straints, bad spectra, or bright foreground stars. We 
make two arrays of square pixels of size 0.025° x 0.025° 
in the (A, r]) sky coordinates covering our analysis region, 
and store the angular selection function calculated using 
the MANGLE routine (Hamilton & Tegmark 2004). At 
the location of each pixel, the routine calculates the sur- 
vey completeness in a spherical polygon formed by the 
adaptive tiling algorithm (Blanton et al. 2003a) used for 
the SDSS spectroscopy. The resulting useful area (with 
nonzero selection function) within the analysis regions is 
1.362 sr. In our study we use only volume- limited sam- 
ples of galaxies, defined by absolute magnitude and red- 
shift limits as shown in Figure 3 of Paper I. In addition 
to the volume-limited samples, labeled Dl to D5, which 
are used to study the environmental effects, we also use 
the galaxy (LS) sample. The LS sample is designed to 
include 'typical' galaxies (—20.0 > Mr > —21.0), and 
covers the full range of comoving distance, from R = 74.6 
to 314.0 h~^M])c, in which our analysis is made. All 
measurements of local density in this paper use the LS 
sample to define the local environment (see section 3). 
The definitions for all samples are summarized in Table 
1. 

We analyze our samples in the survey coordinate 
frame, and the Cartesian coordinates are chosen to be 

X = —RsmX,y = R cos A cos r],z = R cos A sin 77. (1) 

The y-axis is at the center of the survey coordinate 
and radially outward, the z-axis is toward the {rj, A) = 
(90°, 0°) direction, and the x-axis is at the (0°,-90°) 
direction. 

2.3. Physical Parameters of Galaxies 

The physical parameters considered in this study are 
absolute Petrosian magnitude in the r-band Mr, mor- 
phology, '^■^{u — r) color, Petrosian radius, axis ratio, 
concentration index, color gradient in "^'^(.g — i) color, 
velocity dispersion, and equivalent width of the Ha line. 
The r-band absolute magnitude Mr is the AB magnitude 
converted from SDSS magnitudes. To compute colors, 
we use extinction (Schlcgel et al. 1998) and if-corrected 
model magnitudes. The superscript 0.1 means the rest- 
frame magnitude i^-corrected to redshift of 0.1 (Blanton 
et al. 2003b). This makes galaxies at ^ = 0.1 have K- 
correction of — 2.5 log(l -1-0.1), independent of their spec- 
tral energy distributions (SEDs). All of our magnitudes 
and colors follow this convention, and the superscript 
will be subsequently dropped. We also drop the +51og/i 
term in the absolute magnitude. We use the luminosity 
evolution correction of E{z) = 1.6{z — 0.1) (Tegmark et 
al. 2004). 

We employ the g- and i-band atlas image and ba- 
sic photometric parameters measured by the Prince- 
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Fig. 1. — Angular definition of the SDSS sample used for our environmental effect study. Solid lines delineate the boundaries of the 
analysis regions in the survey coordinate plane (A,??). In addition to these boundaries the angular selection function in each spherical 
polygon defined by spectroscopic tiling is fully taJjen into account in our analysis. We name the lower and upper areas regions 1 and 2, 
respectively. 



TABLE 1 
Volume-limited Samples 



Name 


Absolute Magnitude 


Rcdshift 


Distance 


Galaxies(ArB'^) 




Dl 


-18.0 > Mr 


0.025 <z < 0.04374 


74.6 <R< 129.9 


20288 (6256) 


3.41 


D2 


-18.5 > Mr 


0.025 < z < 0.05485 


74.6 < R < 162.6 


32543 (11341) 


3.78 


D3 


-19.0 > Mr 


0.025 < z < 0.06869 


74.6 < R< 203.0 


49571 (19270) 


4.18 


D4 


-19.5 > Mr 


0.025 < 2 < 0.08588 


74.6 < R< 252.9 


74688 (33039) 


4.58 


D5 


-20.0 > Mr 


0.025 < 2 < 0.10713 


74.6 <R< 314.0 


80479 (39333) 


5.56 


LS 


-20.0 > Mr> -21.0 


0.025 < 2 < 0.10713 


74.6 < R < 314.0 


66176 (30715) 


5.94 



^Comoving distance in units of h ^Mpc 

''Number of early type galaxies 

•^Mean separation of galaxies in units of /i~^Mpc 



ton/NYU group ^ to measure the key physical param- 
eters of the 317,533 galaxies in our sample. The g — i 
color gradient is given by the color difference between 
the region with R < 0.5i?pot and the annulus with 
0.5i?pct < R < Rpct, where i?pot is the Petrosian radius. 
Wo use the isophotal position angle, the secing-corrcctcd 
isophotal axis ratio in the i-band, and elliptical annuli 
in all parameter calculations to take into account flat- 
tening or inclination of galaxies. The (inverse) concen- 
tration index, Qn, is defined as R50/R90 where i?5o and 
i?90 are the radii from the center of a galaxy contain- 
ing 50% and 90% of the Petrosian flux in the i-band, 
respectively. The g — i color gradient and concentration 
index are corrected for the effects of seeing using Sersic 
model fits of the g- and i-band images as described in Pa- 
per I. We choose the concentration index as a measure 
of surface brightness profile because the Sersic indices 
measured from SDSS images are quite inaccurate except 
for large bright galaxies with r < 16, due to the rela- 

^ http://photo.astro.princeton.edu 



tively large seeing(~ 1.4") of the SDSS. If the profile fit 
is made by avoiding the central part of galaxy images, 
which is much affected by the seeing, the seeing effects 
on the Sersic index can be reduced but the fit becomes 
more unreliable due to the smaller fitting interval and 
noisier profile. In contrast, by definition Cm = R50/R90 
is insensitive to the surface brightness profile near the 
galaxy center within radius of i?5o. Our seeing-corrected 
concentration index seems accurate because it shows no 
dependence on either redshift or seeing. 

We classify morphological types of galaxies using the 
prescription of Park & Choi (2005). Galaxies are divided 
into early (ellipticals and lenticulars) and late (spirals 
and irregulars) types based on their locations in the u — r 
versus g — i color gradient space and also in the concen- 
tration index space. The resulting morphological clas- 
sification has completeness and reliability reaching 90%, 
as claimed by Park & Choi (2005). When photometry 
is excellent and galaxy images are well-resolved, more 
information from surface brightness fiuctuations can be 
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added for morphology classification. But this is certainly 
not the case near the faint limit (r = 17.6) of the SDSS 
sample we use. For the volume-limited sample D2 we 
perform an additional visual check of the color images of 
galaxies to correct misclassifications by the automated 
scheme for about 20,000 blue early types (those below 
the straight line in Fig. 7 of Paper I) and red late types 
(those with u — r color redder than 2.4). The morpho- 
logical types of 1.9% of galaxies, which are often blended 
or merged objects, are changed by this procedure. 

We use Petrosian radius to represent the physical size 
of galaxies. The Petrosian radii in our analysis, calcu- 
lated by using elliptical annuli, are usually larger than 
those in the DR5 catalog, which adopted circular annuli. 
The angular size is converted to a physical size in units 
of /i~^kpc from the redshift assuming a ACDM universe 
with = 0.73 and f7„ = 0.27. 

The stellar velocity dispersion of SDSS galaxies is mea- 
sured by an automated pipeline called IDLSPEC2D version 
5 (D. J. Schlegel et al. 2007, in preparation). Galaxy 
spectra of SDSS galaxies are obtained by optical fibers 
with angular radius of 1.5". The finite size of the sam- 
pled light smooths the central velocity dispersion profile. 
To correct the central velocity dispersion for this smooth- 
ing effects we have adopted a simple aperture correction 
formula Ucorr = crfiber(8-Rfiber/-Ro)'''°^ (J0rgensen et al. 
1995; Bernardi et al. 2003a), where Rq is the equivalent 

circular effective radius ((6/a)^{v''deV and rdev is the 
seeing-corrected effective angular radius in the i band. 

3. LOCAL DENSITY ESTIMATION 

To find relations between intrinsic physical properties 
of galaxies and their environment we require a well de- 
fined and robust measure of environment. We seek to 
use an environmental parameter that is defined directly 
from observational data, is continuous, and character- 
izes the full range of galaxy environments, from the most 
massive clusters to voids. The local number density of 
galaxies at a given smoothing scale is such an example 
and will be the environmental parameter in our study. 
Local density gradient and shear can be also used to de- 
fine environment, but will be considered in later studies. 
Our choice excludes parameters like the discrete environ- 
mental types such as cluster, group, or field. Assignment 
of galaxies to such environmental types are arbitrary to 
some extent and can be constructed based on continuous 
parameters when desirable. 

Local density estimation depends on the shape and size 
of the smoothing kernel. We adopt the spline kernel with 
adaptive smoothing scale, hg to include a fixed number 
of galaxies within the kernel weighting. This adaptive 
smoothing kernel is often used in smoothed particle hy- 
drodynamics simulations. We use the spline kernel be- 
cause it is centrally weighted, unlike TopHAT or cylin- 
drcal, and has a finite tail, unlike the Gaussian. The 
spline kernel was first used to estimate the local galaxy 
density by Park et al. (1994) who demonstrated that the 
under-dense regions in the Center for Astrophysics (CfA) 
survey lack bright galaxies, indicating an environmental 
dependence of the luminosity function (LF). Details of 
our method to calculate the local density using the spline 
kernel are summarized by C. Park et al. (2007, in prepa- 
ration), who demonstrate that the smooth density field 
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Fig. 2. — Distributions of the local density {top) and the effective 
Gaussian smoothing length estimated at the location of galaxies 
in the D3 sample (bottom). The number of L. galaxies included 
within the spline kernel are 20 {thin histogram) and 200 {thick 
histogram). 

constructed by the spline kernel in redshift space recovers 
the real space density field much better than the cylin- 
drical kernel at low- and intermediate-density regions. 
The local density of galaxies sitting at very high density 
regions and with large peculiar velocities (i.e. those ap- 
pearing at the tips of fingcrs-of-god in redshift space) is 
better meastired by the cylindrical kernel, but the frac- 
tion of such galaxies is very small. 

Because this smoothing kernel is adaptive, it has the 
important feature that it does not oversmooth dense 
regions, thereby blurring the "walls" into the "voids." 
Cosmological simulations clearly show that the transi- 
tion from the interior of large cosmic voids, which have 
nearly constant density, to the denser structures, defined 
by filaments and clusters, is quite abrupt. Maps of the 
observed redshift space distribution of galaxies support 
this picture (see Fig. 3, below). An adaptive smoothing 
kernel, such as the one we employ, better preserves that 
property of the galaxy distribution than a kernel with a 
fixed smoothing length or with an infinite tail. 

On the other hand, the adaptive nature of the smooth- 
ing kernel potentially complicates the interpretation, be- 
catisc the physical smoothing scale systematically varies 
with environment, from a small smoothing scale in clus- 
ters to a very large smoothing scale in voids. However, 
the relationship between environment and smoothing 
scale is nearly monotonic and, unlike fixed-scale meth- 
ods, the adaptive kernel has the advantage that the S/N 
for estimating density is more uniform. 

We define the local density at a position in space as 
the spline kernel-weighted number density of L* galax- 
ies surrounding that position. In every case we use the 
volume-limited LS sample of galaxies with absolute 
magnitudes —21.0 < Mr < —20.0 for local density esti- 
mation. In this way the local density has a clear meaning, 
and it becomes possible to compare the local densities 
of galaxies in different volume-limited samples sitting at 
different radial distances. The number of galaxies re- 
quired within the smoothing volume is set to Ng = 20, 
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which we find is the smallest number yielding good local 
density estimates. Unless otherwise stated, Ng = 20; be- 
low we also examine results with Ng = 200 when study- 
ing larger scales. 

When searching for the spline radius containing the re- 
quired number of galaxies, the smoothing volume can hit 
the sample boundaries or masked regions due to bright 
stars or galaxies. In addition to this trouble, the angular 
selection function varies across the sky even within the 
survey region. We first check if the location in which 
the local density is being measured is within the survey 
area, and then ignore the survey geometry and angular 
selection fimction to search for the nearest 20 galaxies 
in the LS sample in order to find the spline radius and 
thus the local density. We then randomly throw 1000 
test particles into the spline radius around the location, 
and calculate the sum of normalized weights of the test 
particles assigned by the spline kernel, the angular selec- 
tion function, and the radial boundaries. The final local 
density is given by the first estimate divided by the sum 
of weights. Our results use only those local density es- 
timates for which the correction factor is smaller than 
2. 

The size of the spline kernel hs can be related with the 
effective Gaussian smoothing scale using the fact that 
the effective smoothing volumes are nh^ and (27r)^ '^i?G''^ 
for spline and Gaussian smoothings. The result is Rq = 
0.584/is. In the case of galaxies in Sample D3, we find 
the median values and 68% limits of hs are 8.1^|;5 and 

21.0;|° or the effective Ra « 4.7t?;|? and 12.3t^ j when 
Ng = 20 and 200 arc used, respectively. The correspond- 
ing distributions of the local density and the effective 
Gaussian smoothing length are shown in Figure 2. 

4. ENVIRONMENTAL DEPENDENCE OF 
PHYSICAL PARAMETERS 

4.1. Morphology 

The morphology-luminosity relation has been exam- 
ined in Paper I. In this paper we examine this relation as 
a function of local density environment, using our adap- 
tive smoothing to label regions of space. The upper panel 
of Figure 3 shows the distribution of early (red) and late 
(blue) type galaxies in the D3 sample in a wedge defined 
by 26.25 < r] < 36.25 and -55.2 < A < 51.6 in region 
2 (see Fig. 1) of the SDSS. The bottom panel shows 
the same galaxies, but distinguished by their absolute 
magnitudes; bright galaxies with Mr < —20 are red and 
fainter galaxies with Mr > — 20 are blue. 

Figure 4 shows distribution of galaxies in a slab with 
constant thickness of 20/i~^Mpc that cuts through re- 
gion 2. Overlaid are isodensity contours estimated from 
the local densities on a uniform grid through the mid- 
dle of the slab. The upper and lower panels use the 
local density calculated using the nearest 20 and 200 
L* galaxies, respectively. This comparison is made to 
study the dependence of environmental effects on scale 
in section 5. Contours delineate isodensity regions of 
galaxies. The thin lines mark the regions where 
the local density is p/p = 5 or 0.1 for A^^ = 20 and 
p/p = 0.3, 2, or 3 for = 200, respectively. The 
thick curves correspond to the local density equal to 
the mean density of the galaxies, which is measured 
to be p = (5.94/i-iMpc)-3 = 4.77 x 10-^{h-^Mpc)-^. 



We clearly see that early types prefer high-density re- 
gions, while late types are more uniformly distributed. 
Note carefully that the galaxies plotted are brighter than 
Mr = —19.0 while the density field is determined by ~ 
galaxies. Throughout this paper we use the LS sample to 
estimate local density, regardless of which volume- limited 
sample is being examined. 

Figures 5 and 6 present the morphology-density- 
luminosity relation of SDSS galaxies. Figure 5 shows 
the fraction of early type galaxies as a function of local 
density when the absolute magnitude is fixed. We first 
remove late type galaxies in each sample whose sceing- 
corrected i-band isophotal axis ratio is b/a < 0.6 in order 
reduce the internal absorption effect on absolute magni- 
tude (cf. Paper I). We then bin galaxies in the abso- 
lute magnitude versus local density space. The bin size 
in this two-dimensional parameter space is Alogp = 0.8 
and AMr = 0.8. We measure the fraction of early type 
galaxies, the median local density, and the median abso- 
lute magnitude of galaxies belonging to each bin, taking 
into account the fact that we have excluded some frac- 
tion of late type galaxies. When we measure the fraction 
of early type galaxies, we need to know the number of 
early and late type galaxies. Since we are only using late 
types with b/a > 0.6, we should infer the total number 
of late types given the number of those with b/a > 0.6. 
Some late types with b/a < 0.6 must have dropped out 
of the volume- limited sample due to the internal absorp- 
tion. In Paper I we showed that the correction factor can 
be reasonably made using the factor found in the faintest 
volume-limited sample of galaxies with Mr < —17.5. The 
error bars are estimated from 16 random subsets of galax- 
ies. 

We find that the early type fraction is a monotonically 
increasing function of local density at fixed luminosity. 
It is surprising to find that this monotonic dependence 
holds over 3 orders of magnitude in local density, from 
the highest down to the lowest density we explore, and 
in every luminosity bin. The early type fraction falls 
well below 20% at extreme low densities. This late type 
dominance at low density is particularly true for fainter 
galaxies. To test for the effects of the finite bin size in 
log/9 and Mr on the relation we also examine the relation 
using smaller bin sizes. The thin solid line between the 
lines of Mr = —19.9 and —20.4 is the morphology- density 
relation of galaxies in the D5 sample with a median mag- 
nitude of Mr = —20.25 using the bin sizes Alogp = 0.6 
and AMr = 0.5. The amplitude and slope of the relation 
is quite consistent with those for which larger bin sizes 
are used; thus, our results are not affected by the finite 
bin sizes. We also examine the effects of including those 
galaxies whose redshifts arc borrowed from the nearest 
galaxy with redshift within 55", and find that the depen- 
dences of morphology fraction on density and luminosity 
are essentially the same. 

Our results on the morphology-density relation are 
consistent with, and extend to a larger range of den- 
sity and luminosity, the findings of several earlier inves- 
tigators. For example, in studies limited to rich cluster 
environments, Oemler (1974) and Dressier (1980) found 
higher fractions of E and SO types (both classified as 
"early type" here) in the centers of clusters. Subsequent 
investigations have verified and refined the understand- 
ing of this effect. As originally noted by Postman & 
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Fig. 3. — Top: Distribution of early and late morphological types of galaxies in the D3 sample in a wedge defined by survey coordinate 
limits of 26.25 < rj < 36.25 and —55.2 < A < 51.6. Early types are red, and late types are blue. Bottom: Distribution of galaxies with 
absolute magnitudes Mr < —20 (red), and > —20 (blue) in the D3 sample. The lower wedge has the same angular definition as the upper 
one. Abell clusters A2197 and A2199 are located at the lower left corner of the wedges. 



Geller (1984) and by a number of more recent investiga- 
tors (Lewis et al. 2002; Goto et al. 2003; Balogh et al. 
2004a; Weinmann et al. 2006; Quintero et al. 2007) this 
segregation is not limited to rich clusters, but occurs in 
poorer groups as well. Figure 5 reflects this fact in the 
strong dependence of the early type fraction on density 
even at the lowest densities. The dependence of the early 
type fraction on absolute magnitude, even at fixed den- 
sity, has been previously noted (Balogh et al. 2004b; 
Kauffmann et al. 2004; Hogg et al. 2004). As we wiU 
discuss later in more detail, the dependence of the late 
type subclass on density has also been seen as a slight 
dependence of the median blue sequence color on envi- 
ronment (Balogh et al. 2004b; Blanton et al. 2007). The 
main difference between our results here and these previ- 
ous studies is that most of the previous studies use either 
visual classification or use simple classifications accord- 
ing to optical color. Although our classification system 
is related to broadband optical color, using the color gra- 
dients and concentrations yields a cleaner sample mor- 



phologically, bridging the gap between the studies based 
on visual classification and those based on purely color 
classification. Another important difference is the use of 
spline kernel to define local density. The cylindrical ker- 
nel, which is frequently used in past studies (Goto et al. 
2003), is very noisy and cannot distinguish among low 
density regions well while the local density estimated by 
the spline kernel holds a good tight correlation with the 
true density field down to extreme void environments (C. 
Park et al. 2007, in preparation). 

Figure 6 presents the morphology-density-luminosity 
relation in a different way. Here we find that the early 
type fraction is a monotonically increasing function of lu- 
minosity at a given local density. In slightly over-dense 
regions, where the statistics are very good, it can be 
clearly seen that the early type fraction rises abruptly 
at magnitudes brighter than —21.2 (see also Fig. 9 of 
Paper I). The slope of the L-a relation also has a break 
at a similar magnitude (see Fig. 14, below). This may 
imply that many early type galaxies more than about 
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Fig. 4. — Distribution of galaxies in a slab with a constant thickness of 20h~^ Mpc from the 2-axis coordinate of 35h~^ Mpc to 55h~^ 
Mpc passing through region 2. Circles are early types, and crosses are late types. The iso-density contours mark the regions with local 
density of L, galaxies of 5, 1 {thick curves), and 0.1, relative to the mean. In the upper panel, the smoothing volume includes 20 L» 
galaxies, while it is 200 in the lower panel, where the thin lines mark the regions with the local density of p/p = 0.3, 2, or 3. The mean 
density of the L, galaxies is p = (5.94/i"^Mpc)^^ = 4.77 X 10~3(/i~^Mpc)~^. 



1 mag brighter than have origins that are different 
from fainter ones. 

It is remarkable that the early type fraction— absolute 
magnitude relation has similar slope for galaxies with lo- 
cal density ranging over nearly 3 orders of magnitude. 
At intermediate magnitudes, near M^, the early type 
fraction is well described by the morphology-density- 
luminosity formula 

f{E) = [0.33 - 0.074(Af,. + 20.3)] (p/p)^-^^'* (2) 

over the whole density range explored. At magnitudes 
fainter than about Af* + 1 the local density dependence 
monotonically steepens at high density, but becomes 
weaker in under-dense regions (see two bottom curves 
of Fig. 5. A result of this is that in very high density 
regions the early type fraction decreases less steeply at 
fainter magnitudes than it does at lower densities (see the 
top curve of Fig. 6). Balogh et al. (2004b) measured the 
fraction of the red population of galaxies within fixed 



absolute magnitude bins as a function of projected lo- 
cal density E5 adopted by Balogh et al. (2004a). Even 
though galaxies are divided by a different criterion and 
the local density measure is different, the local density 
dependence they found is qualitatively quite similar. We 
emphasize that it is important to exclude highly inclined 
galaxies in this analysis. Otherwise, the number of late 
type galaxies brighter than is significantly underes- 
timated, and the slope of the early type fraction versus 
absolute magnitude relation is systematically biased high 
at bright magnitudes. 

Next we look for dependence of subclasses of late type 
galaxies on environment. Subclass classification of late 
type galaxies can be made in the color-color gradient 
space. We divide the late type galaxies into two sub- 
classes; late late types with u — r < 1.8 and early late 
types with u — r > 1.8. This divides late type galaxies 
roughly into those of Sc type and later and those earlier 
than Sc, respectively (Park & Choi 2005). The bin sizes 
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Flc;. 5. — Fraction of early type galaxies as a function of local 
density at a fixed absolute magnitude (number assigned to each 
line). The magnitude bins used to draw the thick lines are — 18.4± 
0.4, -18.9 ± 0.4, -19.4 ± 0.4, -19.9 ± 0.4, -20.4 ± 0.4, -21.2 ± 0.4, 
and — 22.0±0.4, from bottom to top, respectively. The local density 
bin size is Alogp = 0.8. The thin line is for Mr = -20.25 ± 0.25 
and bin size density Alogp = 0.6. 
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Fig. 6. — Early type fraction as a function of absolute magnitude. 
Labels on solid curves indicate local density, as in Fig. 5. The 
long-dashed, short-dashed, and dotted lines are the morphology- 
luminosity relations obtained from the volume-limited samples Dl, 
D2, and D3, respectively, averaged over all environments. 

are again AM,. = 0.8 and Alog(p/p) = 0.8, and only the 
late types with axis ratio greater than 0.6 are used. In 
Figure 7 we find that the fraction of late spiral types with 
u — r color bluer than 1.8 increases strongly as the lumi- 
nosity decreases. This can be easily understood from the 
color-luminosity relation of late type galaxies (see Fig. 
3a of Paper I). This figure also shows that the late spi- 
ral population increases in low density regions for bright 
spirals, in particular. Even though our sub-spiral class 
division is very crude, it can be seen that the sub-type 
of spiral galaxies does depend on local density at fixed 
luminosity. 

4.2. Luminosity 
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Fig. 7. — Fraction of late types (of late types, not of all galaxies) 
that arc classified as late spiral type galaxies (« — r color bluer than 
1.8), as a function of local density at fixed absolute magnitude (the 
number attached to each line). The bin sizes are AMr = 0.8 and 
Alog(p/p) = 0.8, and points are at the median density and the 
median magnitude of galaxies belonging to e3x:h bin. 

Figure 8 shows the dependence of luminosity of galax- 
ies in the D4 sample on the local density. The upper 
panel is for the early types and the lower is for the late 
types with b/a > 0.6. It is clear that high density regions 
contain very bright galaxies while low density regions do 
not. This indicates that there arc systematic effects of 
local density on the maximum luminosity that can be 
expressed by a change of luminosity function. This en- 
vironmental dependence of the luminosity function was 
found by Park et al. (1994) in the CfA redshift survey 
data, and is now confirmed by the SDSS data over 4 or- 
ders of magnitude of local density. Also evident in Figure 
8 is that this local density dependence of luminosity is 
stronger for the early type galaxies. The dependence is 
weaker but can be still observed for the late types. 

The distribution of galaxies in the vertical direction at 
a given local density in Figure 8 reflects the environment- 
dependent luminosity function which we plot in Figure 
9a. The top points and curve arc the luminosity function 
measured from the D4 sample {Mr < —19.5) and its best- 
fit Schechter function, respectively. Curves lying below 
are the luminosity functions of galaxies located in bins of 
different local density, normalized to the whole volume 
of D4, so that the sum of these curves gives the total 
luminosity function. We exclude late type galaxies with 
6/a < 0.6 in the LF measurements to reduce the bias due 
to the internal absorption, and we weight the remaining 
late types by the inverse of their fraction. The LF at 
the highest density {lower solid line with open circles) 
extends to magnitudes much brighter than those reached 
by the LF in the lowest density environment {dotted line 
with triangles), in agreement our qualitative conclusion 
from Figure 8. To answer the question of whether or not 
the dependence of the LF of galaxies on local density is 
due simply to the dependence of morphology fraction on 
local density or to a change of the shape of the LF. we 
normalize the LFs at the faint end and compare their 
shapes. In each density bin used in Figure 9o we divide 
galaxies into early and late type groups, and recalculate 
their LFs. Figure 96 and 9c show that the shape of the 
LF of each morphological type does change as a function 
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Fig. 8. — Local density dependence of luminosity of early and late 
type galaxies in the D4 sample. Very bright galaxies inhabit only 
high density regions while faint ones populate all densities. The 
result is the local density dependence of the luminosity function. 




0.5 1 5 10 

p/p 

Fig. 10. — Local density dependences of the Schechter function 
parameters M* and a calculated from the D3 sample Filled circles 
are the results from the whole D3 sample (with the inclined late 
types excluded), and open circles and crosses are those from the 
early and late type subsets, respectively. The symbols {squares 
and stars) without connecting lines are those from the D2 sample 
{Mr < -18.5). 
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Fig. 9. — (a) Luminosity functions of all galaxies in the D4 sample 
{upper filled circles) and of those located at different local density. 
The ranges of logp/p are given in the legend at the lower right 
corner. (6) Luminosity functions of early type galaxies at different 
local densities normalized to a value at the faintest magnitude. 
Points and lines have the same meaning as in (a), (c) Same as (6) 
but for late type galaxies 



of local density. In particular, the LF of early types is 
more sensitive to environment than that of late types. 

Therefore, the local density dependence of the LF is due 
both to the change of the early type fraction and that of 
the shape of the LF of each morphological type. 

The visual results of Figure 9 can be quantified by com- 
paring the best-fit Schechter function parameters. We 
divide the galaxies in the D3 sample (Mj. < —19.0) into 
late and early types and into seven different local den- 
sity bins. The solid lines with filled circles in Figure 10 
are the Schechter function parameters M* {upper panel) 
and a {lower panel) when the morphology subsets are 
combined. We use the MINUIT package of the CERN 
program library^ to estimate the parameters by the max- 
imum likelihood method. The estimated Schechter func- 
tion parameters are listed in Table 2. The solid line with 
open circles and the dashed line with crosses in Figure 10 
are those of early and late types, respectively. The lumi- 
nosity of early type galaxies as measured by the charac- 
teristic absolute magnitude, Af», starts to rise steeply for 
local density above the mean density while that of late 
types rises below the mean density, stays more or less 
constant at intermediated density, and then rises again 
in high density regions. Previous work has also found 
that the high-luminosity turnover in the luminosity func- 
tion strongly depends on local density or cluster richness 
(e.g., Lumsdcn et al. 1997; Driver et al. 1998; Marinoni 
et al. 2002; De Propris et al. 2003; Balogh et al. 2004a; 
Hoyle et al. 2005; Popesso et al. 2005; Zandivarez et 
al. 2006). Croton et al. (2005) measured the Schechter 
function parameters for early and late type galaxies in 
the 2dFGRS. Contrary to Croton et al., we find that M* 
for early types is changing slowly below the mean density 
and that luminosity for late types rises at high density 
{p/p > 5). The difference seems to be due to difference 

® http://wwwasdoc.web.cern.ch/wwwasdoc/minuit/minmain.html 
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in morphology classification. They used a morphological 
classifier based on fiber spectra of galaxies. Morphol- 
ogy classifiers based simply on total color or spectrum 
tend to misclassify red spirals into early type, and blue 
ellipticals into late type. The color gradient and con- 
centration constraints in our morphological classification 
strongly reduce such misclassifications. Another impor- 
tant difference in our analysis is the exclusion of inclined 
{b/a < 0.6) late type galaxies. When they are included, 
the number of relatively fainter galaxies is overestimated 
because many of the inclined galaxies are actually bright, 
but appear fainter due to internal absorption. In partic- 
ular, this systematic effect of internal absorption makes 
the a parameter of late types more negative and makes 
the M* of late types appear fainter. The rise of the char- 
acteristic luminosity of late types in high-density regions 
would not be detected unless the bright red spirals, which 
are abundant at high density, are correctly classified and 
the internal absorption effects are eliminated. 

The faint-end slope, a, does not show a monotonic be- 
havior for late type galaxies, consistent with the results of 
Hoyle ct al. (2005) and Croton et al. (2005). It is inter- 
esting that a for early types decreases (steeper faint end 
slope) at local density both below and above the mean. 
At local density near the cosmic mean, the early type 
LF has an exceptionally flat faint end slope: a ~ —0.3. 
Detection of this behavior requires careful morphological 
classification, because faint early type galaxies reside pre- 
dominantly not only in massive halos (which prefer high 
density environment) but also in low mass halos present 
in under-dense regions. The latter tend to be fainter 
bluer early type galaxies that are correctly classified by 
our morphological classifier. If galaxies were classified 
by color or their spectral properties alone, these blue 
early types would have been classified as late types. We 
emphasize again that the major differences between our 
analysis and the previous ones lie in the morphological 
classification and exclusion of highly inclined late type 
galaxies, as well as the new local density estimator. On 
the other hand, one should be cautious about the se- 
lection bias when the abundance of different species of 
galaxies are compared. In the case of the SDSS catalog 
there might be some bias due to the fact that some very 
low surface brightness galaxies are missed. 

The Schechter function parameters are measured from 
our volume-limited samples with magnitude limits of 
-19.0 ± 1.0. In Figure 9 we used the D4 sample (Mr > 
— 19.5) to better see the behavior of the bright tail of the 
LF, and in Figure 10 the D3 sample (A'/,. > —19.0) was 
adopted because this sample has the minimum uncer- 
tainties when both M* and a are concerned. However, 
these samples whose ranges in absolute magnitude are 
only about 1 mag below M,, are not really proper for 
giving accurate estimates of a. The limited range in ab- 
solute magnitude also develops some correlation between 
the estimates of M* and a. This is particularly so at low 
density, where is fainter, and for late types for which 
the faint end slope of the LF is steeper. We compare the 
local density dependences of LF parameters measured 
from all five volume-limited samples with one another to 
see the effects of finite magnitude range. In Figure 10 the 
results from the D2 sample {Mj. < —18.5) are plotted as 
symbols without connecting lines for a comparison. We 
confirm that the behaviors of M* of early and late type 



galaxies as a function of local density were basically the 
same for all samples. It was also true for the a param- 
eter of early types. However, the dependence of a of 
late types on local density showed the systematic trend 
that a became more negative as the faint limit of the 
sample became fainter. In order to measure a of late 
types accurately, one needs to use fainter galaxies, and 
for a fiux-limited sample this makes one rely exclusively 
on the very close faint galaxies. It should be noted that 
the nearby universe is a relatively under-dense part of 
the universe, and that the resulting parameters can be 
also biased due to the local density dependence of the 
LF. Dependence of a on environments is still a matter 
of some debate (e.g., de Propris et al. 1995; Seeker et al. 
1997; Valotto et al. 1997, 2001; Trentham 1998; Andreon 
& Cuillandre 2002; Popesso et al. 2006). 

4.3. Color 

It is well known that galaxy colors are redder in re- 
gions of higher local density (Hogg et al. 2004; Blanton 
et al. 2003c; Tanaka et al. 2004). In this section we 

show this trend is dominated by the dependence of color 
on luminosity and morphology. For this purpose it was 
necessary to examine the dependence of color on local 
density when both morphology and luminosity are fixed. 
We use u — r color in this work because it is a measure 
of the star formation activity of galaxies in the recent 
past. The top panels of Figure 11 show scatter plots of 
u — r color of early {left panels) and late {right panels) 
type galaxies in the D2 sample whose absolute magni- 
tudes are between = —18.5 and —19.3 (the median 
magnitude is —18.9). The thick blue solid lines delineate 
the most probable (mode of the distribution of color at 
fixed absolute magnitude and at fixed local density in- 
terval) color of early types and the median color of late 
types. Late type galaxies with axis ratios b/a < 0.6 
are again excluded. The short-dashed green lines are 
the most probable (early) and median (late) colors of 
galaxies with —19.3 > Mr > —20.1. The correspond- 
ing scatter plots are not shown to avoid confusion. The 
long-dashed magenta line and the red solid line are for 
galaxies in magnitude bins of —20.1 > M,. > —20.9 and 
-20.9 > Mr > -21.7 with median magnitudes of -20.4, 
and —21.1, respectively. 

We first note that there is no break in color as a func- 
tion of local density for either early or late type galaxies. 
Tanaka et al. (2004) reported detection of a break in 
color at a density of logSs ^ 0.4 galaxies ft.y5Mpc^^ (see 
Tanaka et al. for the definition of projected local density 
S5) for galaxies with magnitudes M* -|- 1 < M^ < M* + 2. 
It is clear that this result is simply due to the fact that 
blue late types tend to avoid in very high density en- 
vironments, and it is the morphology bias rather than 
the color bias that causes the effect. This feature in the 
morphology bias appears as a break in the morphology- 
luminosity relation at an absolute magnitude of about 
M* - 1 (Paper I). 

Second, we note that the color of galaxies is a very in- 
sensitive function of local density once the morphology 
and luminosity are fixed. As the local density varies from 
p/p w 15 to about 0.37, the color changes only by 0.03 
for early types and by 0.11 for late types; both get redder 
at higher density. This behavior is insensitive to lumi- 
nosity, which can be seen from the fact that the lines are 
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Fig. 11. — Local density dependence of (a) u — r color, (b) A{g — i) color difference, and (c) cjn (inverse) concentration index. Dots are 
galaxies in the D2 sample with —18.5 > Mr > —19.3. The early type galaxies are on the left, and late types are on the right. Blue solid 
lines are the median curves except for the early types in panels (a) and (b) where the most probable values are drawn. Other colored lines 
show local density dependences at brighter magnitudes of Mr ~ —19.6, —20.4, and —21.1. Uncertainty limits are attached to the relations 
for the faintest and the brightest subsamples. 



nearly parallel to one another. On the other hand, the 
color changes by 0.22 and 0.75 as the magnitude changes 
by 2.2 from Mr = -18.9 to -21.1 for early and late 
types, respectively. The difference between the typical 
colors of early and late types is as large as 1 mag when 
Mr = —18.9. Therefore, the dependence of color on local 
density is dominated by the dependences of color on mor- 
phology and luminosity, which in turn are strongly de- 
pendent on local density, as shown in sections 4.1 and 4.2. 
Figure 11 also shows that blue early types live mainly 
at intermediate- and low-density environments, and that 
some red early type galaxies still exist in the extreme low 
density regions with p/p < 0.1. Tanaka et al. (2004) re- 
ported that the star formation rate indicated by the g — i 
color and the equivalent width of the Ha emission line is 
suppressed in dense environments. It turns out that most 
of the effects they detected are due to the dependence of 
morphology and luminosity on local density. 
To inspect the collective behavior of color of each mor- 



phological type we examine the early and late type se- 
quences in the color-magnitude diagram. The points in 
Figure 12 are early (left panels) and late {right panels) 
type galaxies in the D3 sample located in high (upper 
panels) and low (bottom panels) density environments. 
In the upper panels, solid lines are the most probable 
color for the early types and the median color for the late 
types. In the bottom panels, the most probable and me- 
dian color curves at high density are replotted as dashed 
lines, for comparison with the curves at lower density. We 
find that the shift of the early type sequence in the color- 
magnitude diagram toward bluer color is only 0.025 mag 
when the local density changes from p/p > 10 to < 0.5. 
A larger shift of 0.14 in w — r color is observed for the 
late type sequence. 

The very slight dependence of the location of the red 
sequence on the density field has been noted in the past, 
originally by Sandage & Visvanathan (1978), but more 
recently by others (Balogh et al. 2004b; Hogg et al. 2004; 
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Fig. 12. Dependence of the red and blue sequences in the 
color-magnitude diagram on the local density. Solid curves are the 
most probable (early) and the median (late) relations defining the 
sequences at different environments. The solid curves in the upper 
panels are plotted again in the lower panels as dashed curves for 
comparison. 

Bernardi et al. 2006; Blanton et al. 2007). The stronger 

dependence of the blue sequence colors on environment 
was pointed out by Balogli et al. (2004a). 

4.4. Color Gradient 

The points in Figure 116 are the color gradient of galax- 
ies in the D2 sample as a function of local density. The 
most probable color gradient of early type galaxies is 
slightly negative (the center is redder) and shows no en- 
vironmental dependence at a given absolute magnitude. 
It is nearly independent of luminosity as well. This ex- 
traordinary result is due to the fact that the internal 
color gradient of early type galaxies is almost indepen- 
dent of luminosity even though the color itself does have a 
monotonic dependence on luminosity (cf. Fig. 3a and 36 
of Paper I) . Late type galaxies have more negative color 
gradients when they are brighter (cf. Fig. 36 of Paper 
I), which explains the dependence of the color gradient 
of late types on absolute magnitude in the right panel of 
Figure lib. This may reflect a change in the bulge-to- 
disk ratio with luminosity. Figure 11a shows that fainter 
late type galaxies are bluer in their integrated color than 
bright ones, and Figure 116 indicates that they have bluer 
centers (thus lower color gradients). The color difference 
weakly depends on local density for faint late type galax- 
ies while no dependence is seen for bright galaxies. 

4.5. Concentration Index 

At fixed color or estimated star-formation history, sev- 
eral previous authors found little dependence of con- 
centration on environment (Hashimoto & Oemler 1999; 
Bernardi et al. 2003a; Kauffmann et al. 2004; Blanton et 
al. 2005a; Quintero et al. 2007). We extend these results 
by showing that even when we use a cleaner definition of 
morphology, the radial profiles of galaxies appear to be 
only weakly related to their environment. 

When morphological type and luminosity are fixed, the 
surface brightness profile, represented by the (inverse) 
concentration index, is nearly independent of local den- 
sity as shown in Figure 11c. The median Cin of bright 



early types with absolute magnitudes within 0.8 mag bins 
shows little change at all local densities explored. There- 
fore, the dependence of surface brightness profile of early 
types on local density is solely due to the dependence of 
the profile on luminosity (cf. Fig. 3c of Paper I), which 
monotonically depends on local density, as demonstrated 
in section 4.2. Faint early types seem to be slightly more 
concentrated at low density. It is difficult to see a depen- 
dence of concentration index of late types on local density 
because the dispersion in Cin is large and the effects seem 
weak. 

4.6. Size 

The size (Petrosian radius) of galaxies is an increasing 
function of luminosity as shown in Figure 4a of Paper I. 
Figure 13a shows that physical size is almost indepen- 
dent of local density when luminosity and morphology 
are fixed. We detect a small monotonic dependence of 
galaxy size on local density for both early and late types, 
except for the very bright galaxies. Surprisingly, galax- 
ies of the same luminosity are smaller and therefore have 
higher surface brightness at low density. Galaxies with 
Mr = -19.7 ± 1.2 located at low density {p/p « 0.4) 
are about 8% smaller in size than those at high density 
{p/p « 15), for both early and late types. The weak 
dependence of size on local density should be accepted 
with caution, because including light from neighboring 
galaxies can in principle create a trend in the sizes with 
density in the direction we see. 

Kauffmann et al. (2004) and Blanton et al. (2005a) 
have both investigated the relationship between galaxy 
size and environment (at fixed luminosity and at fixed 
star- formation history or color, respectively), and both 
found only small effects. In the range of luminosities 
we consider here, more luminous than M*, Blanton et al. 
(2005a) find a trend for larger galaxies to be in denser re- 
gions as well (although they found this trend is reversed 
at lower luminosities). Bernardi et al. (2003b) investi- 
gated the dependence of residuals of early type galax- 
ies from the fundamental plane on environment and find 
that the residuals from the fundamental plane in effective 
surface brightness correlate slightly with environment. 
The sense is that at fixed luminosity and velocity disper- 
sion, galaxies are larger in dense regions, similar to our 
results here. 

4.7. Velocity Dispersion of Early Type Galaxies 

In the right panel of Figure 136 we see that the cen- 
tral velocity dispersion of late type galaxies shows no 
dependence on local density. However, we find an in- 
teresting dependence of the L-a relation on environment 
for early type galaxies. The velocity dispersion of fainter 
early types slightly increases toward low-density regions, 
while the opposite is true for brighter early types. This 
might be evidence for an environment-dependent Faber- 
Jackson relation or, more generally, such a dependence 
of the L-(T relation. 

To explore this possibility further we plot early type 
galaxies in high and low-density regions in the abso- 
lute magnitude versus velocity dispersion space. The left 
panels of Figure 14 show the scatter plots of early type 
galaxies in the D4 sample located at different density and 
their most probable relations in the velocity dispersion- 
magnitude diagram. The dashed curves in the lower 
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Fig. 13. — Same as Fig. 11, but for (a) the seeing-corrected elliptical Petrosian radius, (fe) the central velocity dispersion, cr, and (c) the 
equivalent width of the Ha line. All lines are the median relations. 



panel are equal to the solid ones in the upper panel. We 
clearly see that relatively bright early type galaxies in 
higher density regions have higher central velocity dis- 
persions than those in under-dense regions. The right 
panels are for the galaxies in the D2 sample and show 
the relations down to fainter magnitudes. These plots 
indicate that the slope of the L-a relation depends on 
local density and that this trend is more evident at faint 
magnitudes. If we fit a power-law relation L (x a'^ to the 
curves of most probable velocity dispersion, we measure 
7 = 1.8 and 2.2 in regions with local densities of > 10 
and < 0.5, respectively, near Mr = —19.4. Thus, not 
only is the relation between luminosity and a not well 
represented by a single power-law (cf. Fig. Ab of Paper 
I) , the best-fitting such relation depends on local density 
of the galaxies. 

4.8. Star Formation Rate 

The equivalent width of the Ha line is measured from 
the SDSS spectra, and is often used as a measure of re- 



cent star formation activity (or nuclear activity). Be- 
cause SDSS spectra arc obtained from the central 1.5" 
radius region of galaxies, this line width is not necessarily 
a representative measure of star formation activity, par- 
ticularly for spirals with a significant bulge component. 
Nevertheless, the star formation activity as measured 
by the equivalent width of the Ha line, is on average 
stronger for faint late type galaxies (cf. Fig. 4c of Paper 
I). We find that Ha line strength varies little with local 
density once the morphology and luminosity are fixed. 
Note that a strong local density dependence of W{Ha) 
can be obtained if all morphological type and luminos- 
ity are mixed because the fraction of late type and faint 
galaxies increases at lower densities. However, we are 
able to detect a weak correlation between Ha line width 
and local density for late type galaxies; the star forma- 
tion activity signaled by the Ha line strength is relatively 
stronger at low density environment at a fixed luminos- 
ity. For galaxies with Mr = — 18.9±0.4 a linear fit to the 
median relation gives \ogW{Ha) = 1.466 — 0.046 logp/p. 
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Fig. 14. — Central velocity dispersion vs. absolute magnitude 
relations at high (top panels) and low {bottom panels) density cnvi- 
ronnicnts. Solid curves arc the most probable relations (the dashed 
curves in bottom panels are the same as the solid curves in the top 
panels). Straight bars illustrate the slopes of the L-rj relation of 
7 = 4 {left) and 7 = 3 {lower right). 

This dependence is weaker for brighter galaxies. This 

variation of star formation activity with environment is 
consistent with the findings of Rojas et al.(2005), who 
compared star formation activity of galaxies in voids with 
that of galaxies of similar luminosity in denser environ- 
ments. 

Other investigators who have explored the dependence 
of Ha emission on density have found that even at fixed 
galaxy concentration there is a strong relationship be- 
tween the fraction of emission line galaxies and density. 
However, Gomez et al. (2003) and Balogh et al. (2004a) 
found, as we do here, that the distribution of Ha equiv- 
alent widths for emission line galaxies are only weakly 
related to density. These results parallel the similarly 
weak relationship between color and environment for late 
type galaxies. The majority of early type galaxies show 
vanishing star formation activity regardless of luminos- 
ity. We note that there is a group of active early types 
located at median or low density environments. 

4.9. Axis Ratio of Early Type Galaxies 

The points in Figure 15 are the isophotal b/a ratio of 
early type galaxies in the D4 sample with absolute mag- 
nitudes between —19.5 and —20.3 as a function of local 
density. Their median relation is the thick solid line. 
Other lines are for brighter galaxies in the same sample, 
and the corresponding scatter plots are not made, in or- 
der to avoid confusion. Brighter early types are rounder 
than fainter ones (see also Fig. 5 of Paper I), which might 
be explained by the tidal interaction theory (cf. Thuan & 
Gott 1977). Once the luminosity is fixed, the local den- 
sity dependence of the axis ratio of early typos is so weak 
that it is not detected in this plot. The bottom panel of 
Figure 15 compares the histograms of the b/a ratio of 
galaxies at very high and very low density environments. 
No statistically significant difference is observed. 

5. DEPENDENCE ON SMOOTHING SCALE 

So far, we have used local density estimated by the 
spline kernel smoothing in which the kernel size varies to 
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Fig. 15. — Top: Isophotal axis ratio of early type galaxies as a 
function of local density. Galaxies are selected from the volume- 
limited D4 sample {Mr < —19.5). Dots are galaxies with absolute 
magnitudes between —19.5 and —20.3, and the solid line delineates 
the median relation for theses galaxies. The other two lines are for 
brighter galaxies. Bottom: Histogram of the axis ratio of early 
type galaxies in the D4 isample located at high and low density 
environments. Only galaxies with absolute magnitudes between 
— 19.5 and —20.3 are used. 



include Ns — 20 galaxies in the LS sample of galax- 
ies. We now inspect the environmental dependence of 
physical properties of galaxies when the smoothing ker- 
nel includes 200 galaxies. Because we have found that 
most galaxy properties are nearly independent of local 
density once morphology and luminosity are fixed, we 
look for the possible dependence of only the morpholog- 
ical fraction and luminosity function on the much larger 
scale density field. The isodensity contours of the up- 
per and lower plots in Figure 4 compare the two density 
fields p2o/p and p20o/p smoothed over = 20 and 200, 
respectively. 

Figure 16 shows the early type galaxy fraction as a 
function of p200 /p at each fixed luminosity with an abso- 
lute magnitude bin size of 0.8. Surprisingly, the nature of 
the relations is basically the same as that found by using 
a smaller smoothing scale of Ng = 20 (although the rela- 
tion seems steeper for fainter galaxies). We have checked 
that we reach the same conclusion even for Ns = 300. 
Figure 17 demonstrates the dependence of LF of galaxies 
on p20o/p- The Schechter function parameters and a 
depend on the larger scale density p200 in the same way 
that they do on the smaller scale density p2o(see Fig. 
10). The characteristic luminosity strongly increases as 
P200 rises. The type-specific parameters also behave in 
the same way. The only noticeable difference is the a 
parameter of the combined sample at the lowest density, 
but the statistical significance of that single point is low. 
It should be noted that the smoothing scale correspond- 
ing to Ng = 200 is not far from the linear scale (the 
rms density fluctuations is about 0.69 while it is 1.8 for 
Ng = 20). The morphology and luminosity of galaxies 
appear to respond sensitively to the linear density field. 
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Fig. 16. Same as Fig. 5, but the spline kernel-weighted back- 
ground density is estimated from the 200 nearest galaxies at 
each galaxy in the D3 sample. 
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Fig. 17. — Schechter function parameters M* and a as a function 
of the large scale background density calculated by using Ns = 200 
in the spline kernel smoothing. Galaxies in the D3 sample are used. 



Since the density field that has Ns = 200 smoothing is 
correlated with that having A^^ = 20, it is not surprising 
to find similar relations between the local density and 
galaxy properties at both smoothing scales. Therefore, 
it is now very important to identify the scale on which 
the large-scale environmental dependence of morphology 
and luminosity of galaxies originates. To answer this 
question we study the dependence of the morphological 
fraction on very small scale environment by extracting 
samples both with and without close neighbors. We se- 
lect "isolated" galaxies brighter than Mr = —19.5 in the 
D2 sample D2 (M^jim = —18.5). We adopt the isolation 
criteria for a target galaxy with Mr that there can be 
no companion galaxy in the absolute magnitude inter- 
val Mr^i + 1 > Mr > Mr,i - 2, with redshift difference 
within 300 km s"-'^, and with separation across the line 
of sight less than 500/i^^kpc for a candidate galaxy i. 
The solid line in Figure 18 is the early type fraction of 
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Fig. 18. — Early type fractions of isolated galaxies (solid line) 
in the D2 sample brighter than Mr = —19.5. as a function of local 
density {Ns = 20). The long- and short-dashed lines are the early 
type fractions of the galaxies in binary systems or those with one 
or more than one companions, respectively. 



isolated galaxies as a function of p2o/p- Since the mor- 
phological fraction of isolated galaxies with no close com- 
panion still correlates tightly with p2o, we conclude that 
the large-scale environment is definitely one of the factors 
that determine the morphology of galaxies. We reach the 
same conclusion when the isolation criteria are changed 
to 500 km in redshift difference and lh~^ Mpc in 
the separation across the line of sight. In Figure 18 the 
long and short dashed lines are the early type fractions 
of the galaxies with only one neighbor (binary) or those 
with one or more companions (unisolated) , respectively. 
We see that galaxies having close companions are more 
likely to be of early type than isolated galaxies, even at 
the same large-scale density. This demonstrates that the 
morphology of galaxies is determined in response to very 
small scale density as well as to large {Ng = 20) and very 
large {Ng = 200) scale background density fields. 

To narrow down the scale responsible for the origin of 
morphology we push our study further using the densi- 
ties at these three scales in Figures 19 and 20. Figure 19 
shows the relation between the densities estimated with 
Ns = 20 and 200 at each galaxy in the D2 sample. Super- 
posed are contours of constant early type galaxy fractions 
of 0.25, 0.30, 0.35, 0.40, and 0.45 from left to right. We 
note that the contours are basically vertical. The result 
shows that the morphology of galaxies sitting at a fixed 
P2o/p is independent of the larger scale density p20o/p 
and that the dependence of morphology fraction on the 
density with Ng = 200 smoothing is merely because of 
its correlation with the density when Ng ~ 20. There- 
fore, the origin of galaxy morphology must be due to the 
environmental conditions at scales much lower than the 
Ng = 200 scale, which is roughly the 12h~^ Mpc Gaus- 
sian smoothing scale. 

Figure 20 compares P2a/p with the distance to the 
nearest neighbor from a galaxy, rncar; representing the 
small-scale density in this logarithmic plot. For each 
galaxy in Sample D2 (Mr,iim = —18.5) that has absolute 
magnitude limits —20.0 > Mr.i > —22.0 we find the near- 
est companion that lies in the range Mr^i -|- 1.5 > Mr > 
Mr,i — 2 and has radial velocity difference less than 500 
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P20/P 

Fig. 19. — Relation between the density p20 estimated from Ns = 20 spline smoothing and the density p200 smoothing. Galaxies are 
those with Mr < —18.5 in the D2 sample. Superposed on the scatter plot are the contours of constant early type fraction of 0.25, 0.30, 
0.35, 0.40, and 0.45, from left to right. 



km s~"'^, which is a Uttle more generous than that used 
above. If luminosity is proportional to mass, 1.5 mag 
fainter luminosity corresponds to 25% of mass. Our cri- 
terion implies that a companion with mass greater than 
25% of a target galaxy's mass is heavy enough to influ- 
ence the morphology of the galaxy. Superposed are again 
contours of constant early type fractions of 0.32, 0.35, 
0.40, 0.50, and 0.55. It shows interesting dependences of 
morphology fraction on both scales. In general, the con- 
tours have positive slopes, which means the early type 
fraction increases as rnear decreases at fixed p2o/p- At 
fixed Tnear the fraction increases as p2o /p increases, which 
is already shown in Figure 18 for rnoar = 0.5ft.~^Mpc. 
The strongest dependence on Tncar is observed when 
''near ~ 0.2h^^ Mpc. The dependence on rnoar is weaker 
when r„oar is higher or lower. In under-dense regions 
below p2o/p < 0-5 the majority of galaxies are isolated, 
with Tnear > 0.5/i~^Mpc, and the morphology fraction 
is much less dependent on the existence of neighboring 
galaxies. In this under-dense environment galaxies are 
mostly of late type and morphology is sensitive to P20/ P- 
The slope of contours depends on our choice of compan- 



ion criteria and the smoothing size used to calculate the 
morphology fraction. But the general trend described 
above remains robust. A previous study analogous to 
ours is Balogh et al. (2004b) who claimed that the frac- 
tion of emission line galaxies depends on both 3.85 and 
0.77/1^^ Mpc Gaussian smoothing scales. Importantly, 
they gave some evidence that the fraction was mostly 
dependent on local density po.77 where pa.ss is high, but 
depended more on background density ps.ss where pa.ss 
is lower. Their claim is in agreement with our findings 
(compare their Fig. 9 with Fig. 20 of this paper). 

6. DISCUSSION 

6.1. Summary 

We have studied dependences of physical properties of 
galaxies in the Sloan Digital Sky Survey (SDSS) on envi- 
ronment. There are several differences between our anal- 
ysis and previous studies. First, we have used smooth 
local density estimated by the spline kernel enclosing 20 
or 200 L* galaxies to define environment. Because the 
kernel used is centrally weighted, finite, and variable in 
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Fig. 20. — Relation between the distance to the nearest companion and the local density with Ns = 20 smoothing. Dots are early (red) 
and late (blue) type galaxies with —22.0 < Mr < —20.0 in the D2 sample. Contours of constant early type galaxy fraction are superposed 
on the scatter plot. 



TABLE 2 

Local Density-dependent Schechter Function Parameters for the D3 sample 



Morphology 


P/P 




M* - 51ogio/i 


a 




Type 




lO-^h^Mpc-s 






10-3/i'^Mpc-3 


all 


all 


18.47 


-20.22 ±0.03 


-0.81 ±0.03 


16.5 




0.37 


3.59 


-19.85 ±0.07 


-0.70 ±0.10 


2.2 




1.39 


7.36 


-20.11 ±0.04 


-0.74 ±0.05 


5.4 




4.59 


7.18 


-20.21 ±0.04 


-0.69 ±0.04 


5.5 




14.85 


1.50 


-20.61 ±0.08 


-0.99 ±0.07 


1.7 


early 


all 


7.11 


-20.23 ±0.04 


-0.53 ±0.04 


5.3 




0.38 


0.88 


-19.95 ±0.13 


-0.56 ±0.16 


0.6 




1.43 


2.65 


-19.98 ±0.06 


-0.13 ±0.08 


1.7 




4.75 


2.92 


-20.24 ±0.06 


-0.47 ±0.07 


2.2 




15.53 


0.79 


-20.68 ±0.11 


-0.87 ±0.10 


0.9 


late 


all 


12.27 


-20.12 ±0.03 


-0.90 ±0.04 


9.4 




0.37 


2.76 


-19.80 ±0.09 


-0.75 ±0.11 


1.6 




1.39 


4.67 


-20.14 ±0.06 


-0.93 ±0.06 


3.6 




4.59 


4.54 


-20.10 ±0.06 


-0.77 ±0.06 


3.3 




14.37 


0.75 


-20.51 ±0.14 


-1.09 ±0.12 


0.8 



fraction of 0.505 of late type galaxies with b/a > 0.6, is used to infer the total number density of late type galaxies. 
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size, it gives us a wide dynamic range in the estimated 
local density over 4 orders of magnitude. Second, we 
have used an accurate automated morphology classifier 
to group large numbers of galaxies into early and late 
types. In the case of the D2 sample with minor correc- 
tion by visual inspection we expect the completeness and 
reliability to exceed 93%. These high quality morphol- 
ogy sets significantly reduce the biases due to misclas- 
sification. Third, we have used a set of volume-limited 
samples defined by absolute magnitude. Environmental 
dependences are studied for these volume-limited sam- 
ples at fixed morphology and luminosity. In this way we 
have avoided confusion due to mixed dependences that 
occur when flux-limited samples are used. We have also 
attempted to reduce the effects due to internal extinction 
by excluding late type galaxies with axis ratios less than 
b/a = 0.6. Thanks to the large size of the SDSS spec- 
troscopic catalog of galaxies, we arc able to measure the 
environmental dependences to reasonable accuracies even 
after dividing galaxies into many subsamples selected by 
morphology and luminosity. Our major findings are as 
follows. 

1. The fraction of early morphological type galaxies 
is a monotonically increasing function of a (well-defined) 
local density and luminosity. The morphology- density- 
luminosity relation, as measured in this work, should be 
a key constraint on galaxy formation models. It remains 
for us to understand what causes the morphology of a 
galaxy at fixed local density. We find that morphological 
fractions vary strongly with local density, but we find a 
mix of both early and late types at all densities. Thus, 
conditions other than local density must contribute to 
determining the morphology of a galaxy at fixed local 
density. The extra cause might be environmental prop- 
erties other than density. The gravitational potential 
shear field is an example of such an environmental prop- 
erty. Morphology can also depend on the merger history 
of each object, which can be stochastic to some degree. 

2. We find that the environmental dependences of 
various galaxy properties are almost entirely due to the 
dependences of morphology and luminosity on environ- 
ment. When morphology and luminosity are fixed, other 
physical properties, such as color gradient, concentration, 
size, velocity dispersion. Ha emission, and axis ratio, 
are nearly independent of local density. Careful inspec- 
tion reveals weak residual local density dependences for 
several properties. We find that the u — r color of late 
type galaxies is monotonically bluer at low density envi- 
ronments, even at a fixed luminosity. Early types show 
much weaker dependence. We also find that bright early 
type galaxies in high density regions have higher central 
velocity dispersions compared to those in low density re- 
gions, but that the opposite is true for fainter galaxies. 
A consequence is a local density dependence of the L-a 
relation. Because galaxy properties strongly depend not 
only on luminosity but also on morphology, it is clear 
that galaxy properties cannot be determined solely by 
dark halo mass. 

Our results, using a new morphological classifier and 
a robust method of local density estimation, are largely 
consistent with several previous studies attempts (Kaufii'- 
mann et al. 2004; Blanton et al. 2005a; Hashimoto & 
Oemler 1999; Christlein & Zabludoff 2005). In each of 
those papers, the authors found a relationship between 



galaxy color and environment at fixed concentration, a 
measure of galaxy structure associated with morphology 
(Strateva et al. 2001). Conversely, those papers showed 
that at fixed color there is no relationship between con- 
centration and environment. Those results regarding 
concentration do not refer to morphology as we define 
it here, which indeed includes u — r color as one of key 
indicators. In Park & Choi (2005) we have shown that 
morphology can be traced not simply by concentration 
alone, but by a combination of concentration, color, and 
color gradient. Once morphology is defined this way, al- 
most any other indicator of type appears to be nearly 
fixed with environment, even as the fraction of galaxies 
of each morphology changes. Indeed, they imply that be- 
cause the morphological mix changes with local density, 
at fixed concentration one should indeed measure a cor- 
relation between color and density. The converse result, 
that concentration shows only a weak relationship with 
density at fixed color, is a consequence of the fact that 
color is the more powerful type indicator. 

3. We show that morphology of galaxies appears to de- 
pend on local density, from very small scales (< O.l/i^^ 
Mpc) up to the Ns — 20 and 200 smoothing scales cor- 
responding to about 4.7 and 12h^^ Mpc effective Gaus- 
sian smoothing scales. It is shown that the environment 
at scales larger than about 12h~^ Mpc do not indepen- 
dently affect galaxy morphology. This gives us an im- 
portant clue to understanding the origin of galaxy mor- 
phology. The direct cause of galaxy morphology is not 
related with the background density of the linear-regime 
scale at the present epoch. We have also shown that 
morphology of bright galaxies depends on both the lo- 
cal density with Ng — 20 smoothing and the distance to 
the nearest companion galaxy. The dependence on r,icar 
becomes strongest when the neighbor distance is about 
200/i~^ kpc at which the combined effects of strong tidal 
force and long duration of interaction on galaxy prop- 
erty seems to reach the maximum. The importance of 
the smooth density with = 20 rises at larger and 
smaller separations. 

6.2. How Is the Galaxy Morphology Determined? 

In the currently popular scenario of galaxy formation 
and evolution the properties of galaxies depend cither 
only on the formation history of their dark matter ha- 
los or also on nurture processes such as ram pressure 
stripping, strangulation, and harassment. In the cold 
dark matter (CDM) cosmogony, a galaxy-scale dark halo 
forms through a series of mergers of subgalactic objects. 
The over-density at the galaxy scale, which is statistically 
related with the larger scale density, will determine the 
final mass. Other physical properties of the galaxy will 
depend on the final mass, but much variation is expected, 
depending on how the final mass is actually assembled. 
Gott & Thuan (1976) proposed that galaxy morphology 
is determined by the amount of gas left over at maximum 
collapse of the protogalaxy, which in turn depends on the 
ratio of the star formation time scale tsf to the collapse 
time Tcoii- They claimed that ellipticals had relatively 
small TSF / TcoU (if a star formation law was adopted) 
and had finished star formation by the time of maximum 
collapse. In additional to this initial condition and sub- 
galactic-scale evolution, the local environment can also 
directly affect galaxy properties. Ram pressure stripping 
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(Gunn & Gott 1972) is a mechanism that removes gas 
and shuts off star formation when a galaxy orbits a hot 
cluster halo. The morphology of a spiral galaxy can be 
transformed to an early type by this process. High-speed 
encounters of galaxies with other halos typically in clus- 
ters cause impulsive heatings, called harassment (Moore 
ct al. 1996), and can transform spirals to early types. 
Another mechanism, strangulation (Balogh ct al. 2000), 
can also transform morphology through a decline of star 
formation rate due to shut-off of the newly accreted gas 
when a galaxy enters a cluster or group environment and 
loses its hot gas reservoir. There is also observational and 
theoretical evidence that the tidal force field in clusters 
can transform the infalling spirals to early types (Moss 
& Whittle 2000; Gnedin 2003). 

It is important to note that all these nurture mecha- 
nisms proposed so far are basically effective in group or 
cluster environments. In previous studies it was claimed 
that galaxy properties such as star formation activity and 
morphology fraction showed discontinuities at a critical 
density for galaxies somewhat fainter than the galaxy. 
This finding was used to support the proposition that 
properties of the relatively faint galaxies are affected by 
morphology transformation mechanisms working in the 
group or cluster environment (Tanaka et al. 2004; Balogh 
et al. 2004). However, we find that there is no feature in 
the galaxy property versus local density relation even for 
such relatively faint galaxies when luminosity and mor- 
phology arc fixed. The only feature is the sharp decrease 
of the late-type fraction above the critical luminosity of 
about Mr = —21.3 in the morphology versus luminosity 
relation (Paper I). This indicates that there seems to be 
morphology and luminosity-determining processes that 
work at all local densities in a continuous way. 

The strong effect of a close companion on galaxy mor- 
phology as indicated by this study, might imply the exis- 
tence of a morphology transformation mechanism acting 
at a distance without a direct physical contact. What 
is important is that this process seems to be working in 
a wide density environment. Figure 20 indicates that 
this effect is maximal for galaxies having a companion 
at a separation of about 200h~^ kpc. This might be be- 
cause the tidal effects are strong enough to significantly 
enhance the star formation rate, and the duration of in- 
teraction is large as well at this particular separation. 
The morphology transformation process caused by this 
rather secular effects in binary or few-body systems can 
be smoothly connected with the more instantaneous tidal 
interactions with the cluster core or local substructures 
in clusters (Moss & Whittle 2000; Gnedin 2003). 

Galaxies can be either isolated or in binary /multiple 
systems at all environments. For example, in Figure 18 
the fraction of unisolated galaxies is about 20% even 
in the low density region with p2o/p ~ 0.5. A similar 
finding was reported in the existence of red galaxies in 
low density regions (Balogh et al. 2004b). Figure 18 
also shows that the early type fraction is systematically 
higher for galaxies located at higher p2o/p no matter 
whether they are isolated or in systems. The isolated 
early types may have been initially born as early types 
or may be final merger products of initially multiple sys- 



tems. 

It is reasonable to believe that our local background 
density p20 calculated by using neighboring 20 L* galax- 
ies is a good indicator of the ratio tsf/tcoU at the time 
of galaxy formation. The merger rate during the sub- 
sequent evolution and the average rates of all the nur- 
ture processes mentioned above are likely to be increas- 
ing functions of p2o- The fraction of galaxies in bi- 
nary/multiple systems are also an increasing function of 
P20- If these are true and if small tsf/tcou or above evo- 
lution processes yield early type galaxies, we can under- 
stand the overall dependence of morphology on p20- The 
increase of the early type fraction in binary or multiple 
systems can be explained by the tidal disturbance effects. 
It remains to be understood why the dependence on the 
nearest companion distance becomes maximum near the 
particular separation of 200h^^ kpc, which is roughly 
the separation between two typical bright galaxies when 
their dark halos are just touching. 
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